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ABSTRACT 

Exploiting the H-ATLAS Science Demonstration Phase (SDP) survey data, we have determined 
the luminosity functions (LFs) at rest-frame wavelengths of 100 and 250 /im and at several redshifts 
z > 1, for bright sub- mm galaxies with star formation rates (SFR) > lOOAfoyr -1 . We find that 
the evolution of the comoving LF is strong up to z k, 2.5, and slows down at higher redshifts. From 
the LFs and the information on halo masses inferred from clustering analysis, we derived an average 
relation between SFR and halo mass (and its scatter). We also infer that the timescale of the main 
episode of dust-enshrouded star formation in massive halos (M H > 3 x 10 12 M ) amounts to - 7 x 10 8 
yr. Given the SFRs, which are in the range 10 2 — 10 3 Moyr _1 , this timescale implies final stellar 
masses of order of 10 11 — 10 12 M Q . The corresponding stellar mass function matches the observed 
mass function of passively evolving galaxies at z > 1. The comparison of the statistics for sub-mm 
and UV selected galaxies suggests that the dust-free, UV bright phase, is > 10 2 times shorter than the 
sub-mm bright phase, implying that the dust must form soon after the onset of star formation. Using 
a single reference Spectral Energy Distribution (SED; the one of the z ~ 2.3 galaxy SMM J2135-0102), 
our simple physical model is able to reproduce not only the LFs at different redshifts > 1 but also 
the counts at wavelengths ranging from 250 /im to ~ 1 mm. Owing to the steepness of the counts 
and their relatively broad frequency range, this result suggests that the dispersion of sub-mm SEDs 
of z > 1 galaxies around the reference one is rather small. 

Subject headings: galaxies: formation - galaxies: evolution - galaxies: elliptical - galaxies: high redshift 
- submillimeter 
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1. INTRODUCTION 

The star formation history in galaxies is one of the key 
issues we have to understand in order to reconstruct how 
the Universe evolved from small matter perturbations 
at the recombination epoch to the present richness of 
structures. 

Star formation proceeds at a different pace for differ- 
ent galaxies, depending on the physical conditions of the 
available gas. Early Type Galaxies (ETGs) and mas- 
sive bulges of Sa galaxies are composed of relatively old 
stellar populations with mass-weighted ages of > 8 — 9 
Gyr (corresponding to formation redshifts z > 1 — 1.5), 
while the disc components of spiral and irregular galaxies 
are characterized by significantly younger stellar popula- 
tions. For instance, the luminosity-weighted age for most 
of Sb or later-type spirals is < 7 Gyr (cf. Bernardi et al. 
2010, their Fig. 10), corresponding to a formation red- 
shift z < 1 . In general, the old stellar populations feature 
low specific angular momentum as opposed to the larger 
specific angular momentum of the younger ones. 

How can these facts be interpreted in the framework 
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of the hierarchical evolution of the dark matter (DM) 
galaxy halos that have proven to be remarkably success- 
ful in accounting for the power-spectrum (or the spatial 
correlation function) of the large-scale matter distribu- 
tion (e.g. Springel et al. 2006)? A widely held view is 
that the merging of halos is also the principal mechanism 
driving the evolution of the visible part of galaxies (see 
Benson 2010 for a recent review). However, several bod- 
ies of evidence are difficult to reconcile with this scenario. 

• ETGs are characterized by old and homogeneous 
stellar populations. Correlations tight enough to 
allow little room for random processes such as a se- 
quence of mergers (apart from small mass additions 
through minor mergers at late epochs, see Kaviraj 
et al. 2008), and sensitivity to the environment 
(color—luminosity; fundamental plane relations; 
dynamical mass— luminosity) have been known for 
a long time and have been recently confirmed with 
very large samples, and shown to persist up to 
substantial redshifts (Renzini 2006; Clemens et al. 
2009; Thomas et al. 2010; Rogers et al. 2010; Pee- 
bles & Nusser 2010, and references therein). More 
recently, a remarkably tight luminosity— size cor- 
relation has been reported (Nair et al. 2010). In 
addition, ETGs were found to host supermassive 
black holes whose mass is proportional to the bulge 
and to the halo mass of the host galaxy (see Magor- 
rian et al. 1998; also Ferrarese & Ford 2005 for a 
review). All that indicates that the formation and 
evolution of ETGs is almost independent of envi- 
ronment, and driven mainly by self- regulation pro- 
cesses and intrinsic galaxy properties such as mass. 

• There are rather tight, albeit not inescapable, 
observational constraints on the star-formation 
timescale in the most massive ETGs. An upper 
limit comes from the observed a-enhancement or, 
more properly, iron under-abundance compared to 
a elements. Depending on the slope of the assumed 
initial mass function (IMF), the observed a/Fe el- 
ement ratios requires star formation timescales < 
10 9 yr (e.g. Matteucci 1994; Thomas et al. 1999). 
But in merger-driven galaxy formation models star 
formation in ellipticals typically does not truncate 
after 1 Gyr (Thomas & Kauffmann 1999; however, 
see Arrigoni et al. 2010, Khochfar & Silk 2010). 
// a standard IMF is assumed a lower limit comes 
from (sub-)mm counts, implying that several per- 
cent of massive galaxies are forming stars at rates 
of thousands M© yr -1 at z ~ 2 — 3 (see Chap- 
man et al. 2003, 2005). This requires that this 
star- formation rate (SFR) is sustained for > 0.5 
Gyr, much longer than the timescale of a merger- 
induced starburst, which is of order of the dynam- 
ical time (~ 0.1 Gyr for the massive ETGs of in- 
terest here; see, e.g., Benson 2010; Hopkins 2011). 
In other words, a single starburst episode is too 
short to account for the space density of (sub- 
)mm bright galaxies as well as for their present- 
day stellar masses. And indeed models envisage 
a sequence of mergers, and associated starbursts, 
throughout the galaxy lifetime, i.e. over a time 
span much longer than the upper limit set by the 



a-enhancement. The problem of accounting for the 
counts of sub-mm galaxies may be eased assuming 
a top-heavy IMF (Baugh et al. 2005). Indeed, re- 
cent observational evidences (Gunawardhana et al. 
2011; Dunne et al. 2011) indicate that highly star 
forming galaxies have IMFs dN/dm oc mT x with 
flatter high-mass power-law slopes x than galaxies 
with low star formation rates. Gunawardhana et al. 
(2011), using a sample of galaxies from the Galaxy 
And Mass Assembly (GAMA) survey (Driver et al. 
2009) covering the redshift range < z < 0.35, find 
a dependence of x on the SFR that, extrapolated 
to an SFR = 1000 M yr" 1 , would give x » 1.5, 
substantially flatter than the Salpeter (1955) slope 
(xs = 2.35). However, it is not clear that these re- 
sults apply to the high-z proto-spheroidal galaxies 
considered in this paper, since the IMF may depend 
on other parameters, such as age and metallicity. 
Moreover, the most recent study of the evolution of 
galaxies in the far-infrared/sub- mm based on star- 
bursts triggered by mergers (Lacey et al. 2010) 
resorts to an even flatter high- mass IMF (x = 1) 
and has still serious problems with reproducing the 
Herschel counts (see §[7|). 

• Integral-field near-IR spectroscopy of galaxies with 
less extreme SFRs (of few 10 2 M© yr -1 ) at z ~ 2, 
that appear to be very productive star formers in 
the universe (Dekel et al. 2009), has shown that 
in many cases they have ordered, rotating veloc- 
ity fields with no kinematic evidence for ongoing 
merging (Genzel et al. 2006; Forster-Schreiber et 
al. 2009). Still, they harbor several starforming 
clumps: a complex morphology is not necessar- 
ily a symptom of merging. These galaxies show 
tight SFR— mass correlations, with small disper- 
sions (Daddi et al. 2007; Pannella et al. 2009; 
Dunne et al. 2009; Rodighiero et al. 2010; Maras- 
ton et al. 2010). This is not easily reconciled with a 
scenario in which star formation proceeds through 
a series of short starbursts interleaved by long peri- 
ods of reduced activity and these galaxies have been 
caught in a special, starburst moment of their exis- 
tence. The data are more easily accounted for if the 
high SFRs are sustained for some 1 — 2 Gyr, much 
longer that a dynamical time typical of starbursts. 
Although the duration of the star-formation phase 
for these objects is longer than that of the more 
extreme objects considered in the previous bullet 
by a factor of 2 to 3, their final stellar mass is sev- 
eral times lower because the SFR is about an order 
of magnitude smaller. For these galaxies a weak 
a-enhancement is predicted, consistent with obser- 
vations. 

• A comparison of the stellar mass functions at z > 
1.5 with the local one shows that little additional 
growth can have occurred for z < 1.5 through mi- 
nor mergers (Mancone et al. 2010; Fan et al. 2010; 
Kaviraj et al. 2009). 

It is clear from the above that the reconstruction of the 
star formation history through cosmic time for galaxies 
of different masses provides a critical test for galaxy for- 
mation and evolution. Since the star formation occurs 
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within dusty molecular clouds and is deeply obscured at 
ultraviolet and optical wavelengths, data at far-IR/ (sub- 
)mm wavelengths, where the absorbed radiation is re- 
emitted, are essential to provide a complete picture of 
it. 

In this paper we focus on high redshift (z > 1) galax- 
ies. We exploit the far-IR/sub-mm data collected by 
the Herschel Space Observatory (Pilbratt et al. 2010) 
during the Science Demonstration Phase (SDP) of the 
Herschel Astrophysical Terahertz Large Area Survey (H- 
ATLAS, Eales et al. 2010a). The H- ATLAS is an open- 
time key program that will survey ~ 550 deg 2 with PACS 
(Poglitsch et al. 2010) and SPIRE (Griffin et al. 2010) 
in five bands, from 100 to 500 /im. 

The observed SDP field covers an area of « 3.8 x 
3.8 deg 2 centered on (a, S) « (09 h 05 m , +0°30'). Com- 
plete descriptions on reduction of PACS and SPIRE data 
are given in Ibar et al. (2010) and Pascale et al. (2011), 
respectively. Source extraction and flux density estima- 
tion are described in Rigby et al. (2011). The 5cr detec- 
tion limits, including confusion noise, are 33.5, 37.7, and 
44.0 mJy/beam in the SPIRE bands at 250, 350, and 
500 fj,m, respectively; in the PACS bands they are 132 
mJy and 121 mJy at 100 and 160 /im, respectively. 

The plan of the paper is the following. The selection 
of the sample is described in § [2] In § [3] we discuss the 
far-IR Spectral Energy Distribution (SED) of high-z star 
forming galaxies, a fundamental ingredient for our pho- 
tometric redshift estimates, presented in § [U In § [5] wc 
estimate the galaxy luminosity functions (LFs) at differ- 
ent redshift in the range 1 — 4. In § [5] we discuss some 
clues on star-formation timescales for massive galaxies. 
In §[7]modcl predictions for source counts from 250 (im to 
2 mm are compared with the data. Our main conclusions 
are summarized in § [SJ 

Throughout the work we adopt a standard, flat ACDM 
cosmology (see Komatsu et al. 2011) with matter den- 
sity parameter VLm = 0.27 and Hubble constant Hq = 
70kms~ 1 Mpc" 1 . We adopt a Bardeen et al. (1986) 
cold dark matter power spectrum with primordial index 
n s = 1 and cosmic mass variance as = 0.81. Stellar 
masses and luminosities of galaxies are evaluated assum- 
ing the Chabrier's (2003) initial mass function (IMF); 
these can be converted to a Salpeter (1955) IMF on mul- 
tiplying by a factor w 1.6. 

2. SAMPLE SELECTION 

The H-ATLAS sources comprise both a low-z galaxy 
population, identified through matching to the Sloan 
Digital Sky Survey (SDSS York et al. 2000) data (Smith 
et al. 2011a), and a high-z population (median redshift 
~ 2) identified through their far-IR colours (Amblard 
et al. 2010). Low-z galaxies are generally normal/star- 
forming late-type galaxies with moderate opacity (Smith 
et al. 2011b; Dunne et al. 2011). Through analyses 
of clustering these two populations are found to be very 
different; the low-z population (z < 0.3) clusters like 
star-forming blue galaxies (Guo et al. 2011; van Kam- 
pen et al. 2011; Maddox et al. 2010), while the high-z 
population clusters much more strongly, suggesting that 
the high-z sources reside in more massive halos (Maddox 
et al. 2010). 

In this work we investigate the evolution with cosmic 
time of high redshift galaxies with intense star-formation 



activity, interpreted as massive proto-spheroidal galax- 
ies in the process of forming most of their stellar mass 
(Granato et al. 2004; Lapi et al. 2006). Since these 
objects are observed to be in passive evolution at z < 1- 
1.5, we confine ourselves to z > 1. At these rcdshifts 
sources above the H-ATLAS detection limits have dust 
luminosities > 10 12 L Q and star-formation rates (SFRs) 
> 100 Mq yr" 1 . They are therefore ultra-luminous in- 
frared galaxies (ULIRGs). Their dust heating mostly 
comes from young massive stars within molecular clouds, 
implying, on one side, that their far-IR SEDs are gener- 
ally (albeit not always, e.g. Hwang et al. 2010) warmer 
than those of low-z dusty galaxies which have higher 
contributions from cooler interstellar dust heated by old 
stars, and, on the other side, that their optical emission 
is strongly attenuated. The last point means not only 
that spectroscopic redshifts are available just for a tiny 
fraction of sources, but also that we do not have at our 
disposal the multi- frequency optical/near- IR photome- 
try that allowed photometric redshift estimates at lower 
z (Dye et al. 2010; Vaccari et al. 2010; Eales et al. 
2010b). 

In fact, for most H-ATLAS z > 1 galaxies, the only 
available data is the Herschel photometry, primarily in 
the three SPIRE bands (250, 350, and 500 //m), plus 
mostly upper limits in the PACS 100 and 160 fim bands. 
A key issue is then whether the redshift estimates that 
can be obtained from such data are sufficient to obtain 
meaningful estimates of the LFs at least over a limited 
redshift range. At first sight one would be inclined to an- 
swer 'no', but a closer investigation can suggest a more 
optimistic conclusion. 

The dust re-radiation in starburst galaxies is expected 
to come from at least three astrophysical settings (e.g. 
Silva et al. 1998): molecular clouds, diffuse low-density 
clouds (cirrus), and circum- nuclear regions, heated by 
Active Galactic Nuclei (AGNs). The AGN dust emission 
peaks in the mid-IR (Granato & Danese 1994; Andreani 
et al. 2010; Lutz et al. 2010; Hatziminaoglou et al. 
2010) and can be safely ignored in the SPIRE wavelength 
range. Molecular clouds are the preferential sites of star 
formation implying that they are endowed with intense 
radiation fields and relatively warm dust temperatures. 
The cirrus component is exposed to the less intense gen- 
eral radiation field due to older stellar populations that 
have come out from their native molecular clouds and is 
therefore characterized by lower dust temperatures. 

In the nearby Universe, molecular clouds and cir- 
rus give comparable contributions to the far-IR emis- 
sion from 'normal' late- type galaxies, with relatively low 
SFRs (Rowan- Robinson et al. 2005). The colder cirrus 
contribution is especially important in less optically ob- 
scured IR galaxies (Hwang et al. 2010) while the warmer 
molecular cloud (sometimes referred to as 'starburst') 
component becomes increasingly important for higher 
and higher SFRs (Rowan- Robinson et al. 2010). This 
argument also highlights a possible degeneracy: a 'cold' 
observed SED may be associated either to a low-z cirrus 
dominated galaxy or to a redshiftcd warm galaxy. If the 
redshift is estimated using a warm SED, cold low-z galax- 
ies would be erroneously assigned high redshifts. This 
problem can be overcome, however, because cold, low-z 
galaxies are only moderately obscured by dust (the cirrus 
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optical depth cannot be very large) , and are therefore rel- 
atively bright in the optical bands. This is illustrated by 
Fig. [TJ which shows that the SEDs of optically identified 
z < 0.5 SDP galaxies studied by Smith et al. (2011b) 
imply r-band magnitudes brighter than the SDSS DR7 
limit (r = 22.4) at all redshifts even if their 250 /jm flux 
density is at the detection limit, while z > 1 ULIRGs, 
with SED like that of Arp220, (and even more, younger 
high-z galaxies with SED like that of SMM J2135-0102, 
'The Cosmic Eyelash'; Ivison et al. 2010a, Swinbank et 
al. 2010) are fainter than this magnitude limit for the 
same 250 /zm flux density 1 . 

Therefore, we may weed out cold low-z galaxies by 
dropping SDP galaxies with SDSS counterparts (Smith 
et al. 2011a), except those with optical (spectroscopic 
or photometric) redshifts in the range of interest here 
(z > 1). This operation, however, has potential draw- 
backs. First, it leaves in low-z ULIRGs. This is not a 
big problem, since these objects have warm SEDs and 
therefore, as discussed below, their redshifts can be es- 
timated with sufficient accuracy. Second, the reliability 
of SDSS counterparts can never be 100%. For example, 
strongly lensed galaxies generally have an apparently re- 
liable counterpart which is most likely the foreground 
lens. These objects can however be recovered since most 
frequently the lenses are ellipticals, whose optical colours 
are incompatible with a large dust emission. The fraction 
of false identifications among the optical counterparts to 
H-ATLAS SDP sources with reliability R > 0.8 is esti- 
mated to be w 5.8% (Smith et al. 2011a). Although 
this fraction is reassuringly small, we need to keep in 
mind that some truly high-z sources can be missed by 
our procedure. 

Most importantly, not all the true r < 22.4 counter- 
parts can be identified with > 80% confidence, due to 
incompleteness of the SDSS catalog, positional uncer- 
tainties, close secondaries, and the random probability 
of finding a background source within that search ra- 
dius (Dunne et al. 2011). According to Smith et al. 
(2011a), about 60% of the 6621 with 250 /im flux density 

> 32mJy have counterparts brighter than r = 22.4 mag 
in the SDSS (and are therefore, with few exceptions, at 
z < 1). Of these, 2423 could be identified with a reliabil- 
ity R > 0.8, implying that another ss 1550 sources, i.e. 
about 23% of the total sample are really at z < 1 but are 
missing a reliable identification. From the redshift dis- 
tribution of reliable identifications (Fig. 6 of Smith et al. 
2011a) we estimate that about 20% of these sources are 
at z > 0.5. If the same proportion applies to unidentified 
sources, the identification incompleteness of z < 0.5 H- 
ATLAS SDP galaxies is « 19%, that we conservatively 
round to 20%. The effect of this incompleteness on our 
luminosity function estimates is discussed in §[S] 

Starting from the catalog of SDP sources by Rigby et 
al. (2011), which contains 6876 sources, we drop the 
galaxies for which Smith et al. (2011a) have identified 
reliable counterparts. We further require detection at 

> 3<T at 350 /im. In this way we get a sample (that will 

1 Note however that a galaxy with the real Arp220 luminosity 
would be brighter than r = 22.4 up to z ss 0.6. At z = 0.5 it 
would have r ~ 21.5 mag. At the same redshift, the true, non- 
demagnified SMM J2135-0102, would have r ks 19.6 mag. On the 
other hand, high-2; moderately obscured galaxies can be detected 
by the H-ATLAS survey only if they have very high stellar masses. 



be taken as our reference sample) defined by the following 
criteria: i) S^so^m — 35 mJy; ii) no optical identification 
with R > 0.8 (Smith et al. 2011a); iii) detection at 
> 3cr at 350 /xm. The resulting sample is made of 3469 
sources. The redshift estimates presented in §2] indicate 
that some of these (376) are z < 1 ULIRGS with SEDs 
akin to our templates; they will be excluded from the 
subsequent analysis. 

Most (2763, i.e., « 80%) galaxies satisfying the above 
criteria are detected at > 4tr at 350 /xm, and this obvi- 
ously helps with the photometric redshift estimates. This 
is because low-z galaxies have low S^o^m /^so^m ratios 
and are therefore lost when we raise the 350 /xm flux limit 
(see also Amblard et al. 2010). 

3. SEDS OF HIGH-Z SUB-MM GALAXIES 

For the objects of interest here, with SFR > 100 Mq 
yr , only the 'warm' (starburst) component is relevant. 
It is important to take into account, however, that its 
SED is much broader than a single temperature grey- 
body (see, e.g., Silva et al. 1998). Over the 50 — 500 /xm 
range (in the rest-frame) such a SED can be approxi- 
mated, to better than 10 — 20%, by a sum of two grey- 
bodies, each described by 



with temperatures rs 30 K and « 60 K, and dust 
emissivity indices j3 = 1.7 and 2, respectively (see also 
Dunne & Eales 2001). A fit to the observed SEDs of 
standard template starburst galaxies with quite differ- 
ent SFRs (M82, Arp220), and of the z sa 2.3 strongly 
lensed galaxy SMM J2135-0102, confirms the validity of 
this approximation and shows that the relative normal- 
ization of the two grey bodies varies by factors of sev- 
eral. Mid/far-IR data on starburst galaxies emphasize 
the warmer component, while sub-mm data emphasize 
the colder one; if a single temperature dust is used, data 
in different wavelength ranges may yield substantially 
different temperatures. If the source redshift is in the 
range 1 < z < 3.5, PACS and SPIRE wavelengths ex- 
tend across the dust emission peak (which is typically at 
a rest-frame wavelength A w 90 — 100 /im). This may 
allow reasonably accurate redshift estimates using only 
Herschel data (Negrello et al. 2010). However, only a mi- 
nor fraction of the H-ATLAS SDP galaxies have at least 
one PACS detection, usually at 160 /xm, and these are 
mostly at low redshifts; we must also beware of the flux 
boosting by confusion, increasing at longer wavelengths 
(see §S]). 

Arp220 and SMM J2135-0102 are of particular inter- 
est because their SEDs are well determined and have 
SFRs quite typical of the galaxies considered here. We 
have modeled their spectral energy distributions from 
extreme-UV to radio frequencies through the spectropho- 
tometric code GRASIL 2 , which includes a sophisticated 
treatment of dust reprocessing (Silva et al. 1998; Schurer 
et al. 2009; Silva et al. 2011). The results are il- 
lustrated in Fig. [2j it must be stressed that the dom- 
inant contribution for wavelengths A > 30 /im is pro- 

2 For information on the code see 
http : //adlibitum. oat . ts . astro . it/silva/def ault .html, and 
for a web-based version see http : //galsynth . oapd . inaf . it 
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vided by molecular clouds (more details in Bressan et 
al. 2011, in prep.)- Assuming a Chabrier (2003) IMF 
we obtain a SFR of ss 14OM yr _:L for Arp220 and of 
» 27OM yr- 1 for SMM J2135-0102 (a Salpeter IMF 
would yield w 224M yr _1 and « 43OM yr _1 , respec- 
tively). Then assuming the usual linear scaling between 
the SFR and the continuum far-IR luminosity Lfir in- 
tegrated between 8 /xm and 1000 pm, i.e., 



SFR = k x 10- 44 LfiR 



M yr- 1 



ergs" 



(2) 



we find /cArp = 2.5 and &smm = ^-4. The difference 
between the two coefficients is well within the expected 
range. Kennicutt (1998) pointed out that variations of 
?» 30% can be due to differences in the star-formation 
history, implying different effective ages of the stellar 
populations. In Fig. [5] we also sketch the far-IR SED 
of G15.141 (H-ATLAS J142413.9+022304), a strongly 
lensed sub- mm galaxy at z ~ 4.24 with estimated SFR of 
several hundreds M yr" 1 (Cox et al. 2011), modeled as 
the sum of two grey-bodies with T\ = 32 K, T2 = 60 K, 
(3 = 2, and a ratio of 0.02 between the coefficients of 
the warm and of the cold components. We choose the 
SEDs of these 3 galaxies (Arp220, SMM J2135-01012, 
and G15.141) as our references to quantify the effect of 
different choices for the SED on estimates of the LFs. 
SMM J2135-01012 has the smallest S^m/Sioo^m ratio 
and G15.141 exhibits the steepest decrease at long wave- 
lengths. 

4. ESTIMATING THE REDSHIFTS OF SUB-MM GALAXIES 

A major source of concern is the boosting of fluxes in 
the Rigby et al. (2011) catalogue due to confusion by 
faint sources. Since the effect is larger at longer wave- 
lengths, because of the poorer angular resolution, it tends 
to bias high the redshift estimates based on Herschel pho- 
tometry. According to the simulations carried out by 
Rigby et al. (2011), 56.5% of sources detected at > 5a 
at 500/im show a flux boosting by a factor > 1.5 and 
27.3% by a factor > 2, and the effect becomes negligible 
for > lOcr sources. 

To investigate the quantitative impact of the flux 
boosting we estimated the photometric redshifts of 39 
H-ATLAS galaxies at z > 0.5 for which spectroscopic 
redshifts are available. The results are presented in Fig. [3] 
for our 3 reference SEDs. The average SED of low- 
z H-ATLAS SDP galaxies, determined by Smith et al. 
(2011b), was also used, for comparison. There is no in- 
dication that photometric redshifts of the high-z sources 
are systematically overestimated when we use the SED 
of SMM J2135-0102 as a template. The median value of 
Az/(1 + z) = (zphot - 2 spC c)/(l + z sp cc) is 0.01 with a 
dispersion of 0.21. The situation is only slightly worse 
in the case of Arp220: the median value of Az/(1 + z) 
is 0.06 with a dispersion of 0.27. The median offset be- 
tween photometric and spectroscopic redshifts increases 
to 0.18, with a dispersion of 0.28, if we use the cooler SED 
of G15.141. In the redshift range (z > 1) we are most 
interested in there are 24 objects. The median values of 
Az/(l + z) are -0.002 (dispersion 0.12) for SMM J2135- 
0102, 0.07 (dispersion 0.16) for Arp220, 0.16 (dispersion 
0.13) for G15.141. There is no statistically significant 
difference between the 6 strongly lensed objects, having 



> lOcr detections at 500/im, and the other 33 sources, 
whose flux densities are representative of those in our 
sample. The median values of Az/(1 + z) and the dis- 
persions rapidly increase as we go down in redshift, as 
expected since our templates do not match those of low- 
z galaxies whose SEDs include large cirrus contributions. 
For the 1096 H-ATLAS SDP sources with z < 0.5 and 
spectroscopic redshift in the Smith et al. (2011a) catalog 
we find median Az/(1 +z) of 0.38 (dispersion 0.45), 0.42 
(dispersion 0.54), and 0.63 (dispersion 0.55) for SMM 
J2135-0102, Arp220, and G15.141, respectively. Not sur- 
prisingly, the redshift estimates based on the mean low- 
z SED of Smith et al. (2011a) go in the opposite di- 
rection: the redshifts are systematically underestimated. 
The mean offset is small (Az/(l + z) = -0.08) at z < 0.5 
and increases (in absolute value) to —0.21 for z > 0.5, 
and to —0.34 for z > 1. The dispersions are 0.31, 0.29, 
and 0.23 for z < 0.5, z > 0.5, and z > 1, respectively. 
This confirms that the low-z, optically bright, H-ATLAS 
galaxy population has far-IR properties different from 
those of high-z galaxies. Hence, the SEDs determined at 
low-z are not applicable at high-z. 

This test suggests that, at least for sources at z > 1, er- 
rors on photometric redshift estimates are more related 
to the choice the SED template than to the signal-to- 
noise ratio at 500/im, and hence to flux boosting. This 
is mostly due to the fact that the relative errors on flux 
densities are larger at longer wavelengths so that the data 
points more liable to flux boosting weight less in the min- 
imum x 2 fit t° the template SED. We do see in most 
cases that the best fit SED is below the 500/im data, as 
expected if the latter is overestimated. We have checked 
this by lowering the 500/im flux densities first by 20%, 
and this left the derived redshift distribution almost un- 
changed, and then by a factor of 2. In the latter case, we 
got a median value of (zdeboost — z )/(l + z ) ~ —0.1, with 
only 5% of cases at < —0.2 (here z is the redshift esti- 
mated using the fluxes tabulated by Rigby et al. 2011). 

A second concern is the effect of the SED variety of 
active star forming galaxies. To quantify the correspond- 
ing errors in the redshift estimates, we generated a cat- 
alogue of 9 x 10 3 galaxies with a distribution of flux 
densities reflecting that of our sample and redshifts ran- 
domly selected in the range 1 < z < 3.5. Each redshift 
was randomly assigned to one of the 19 well observa- 
tionally determined SEDs of local (ultra-) luminous IR 
galaxies with star formation rates > 20 M yr -1 and 
with contributions from an active nucleus to the far-IR 
flux /agn < 10% studied by Vega et al. (2008). To the 
250, 350, and 500/im flux densities we associated errors 
extracted randomly from the distribution of the errors for 
real observations. We have then estimated the redshifts 
of simulated galaxies, based on flux densities at SPIRE 
wavelengths, using each of our three template SEDs. 

In w 95% of the cases we have |z p hot — 2true|/ (1 + z) < 
0.2. Note that this is likely an upper limit since the 
previous test and the analysis of multi-frequency source 
counts (see §[7]) strongly indicate that the SEDs of high- 
z galaxies with intense star-formation are more uniform, 
and closer to that of SMM J2135-01012, than those of 
star-forming low-z galaxies. This can be expected since 
these galaxies have generally lower star-formation rates 
(of order of tens M /yr, to be compared with > 100 M 
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yr _1 for high- z galaxies) and higher contributions to dust 
heating from old stellar populations. As the SFR in- 
creases, the ratio between the contributions to the far- 
IR/SMM SED from the molecular cloud component (as- 
sociated to star-formation) and from the cirrus compo- 
nent (heated by the general interstellar radiation field) 
increases causing a systematic variation of the SED shape 
with the specific star formation rate, as indeed observed 
(da Cunha et al. 2008; Smith et al. 2011b). This trend 
stops when the former component dominates, as in the 
case of sources considered in this paper, and the SEDs 
become more uniform. 

We have estimated the redshifts 3 of objects in the sam- 
ples defined in §[2] Sources undetected at > 3 a at 500 //m 
were attributed a 3a upper limit of 27 mJy. The calcu- 
lations were repeated using 5er upper limits (45 mJy): 
the derived redshift distribution did not change appre- 
ciably. The results are similar, but somewhat more sen- 
sitive to the effect of boosting (based on the tests de- 
scribed above), if we use the catalogued flux densities 
and errors also for sources with S^oOfim < 3<r. The red- 
shift estimate is the result of a minimum x^-fit of each 
of the three templates to the SPIRE and PACS data 
(including 3er upper limits 4 ). Fig. @] shows that, after 
correcting for the offsets highlighted by Fig. [31 the de- 
rived redshift distributions are only moderately affected 
by the choice of the template SED; they have broad max- 
ima in the range 1.5 < z < 2.5 and a tail extending up 
to z k. 3.5, consistent with earlier estimates for BLAST 
(Ivison et al. 2010b) and Herschel (Amblard et al. 2010; 
Eales et al. 2010) samples, when the different selec- 
tions, and in particular the fact that we have dropped 
galaxies with R > 0.8 optical counterparts, are taken 
into account. Sources with estimated z < 1 will be ex- 
cluded from the subsequent analysis. With our preferred 
SED, that of SMM J2135-01012, our reference sample 
(<S250j*m > 35 mJy, S 350l j, m > 3cr, and no R > 0.8 optical 
identifications) contains 3093 galaxies (i.e. w 45% of the 
6876 H- ATLAS SDP galaxies) at z > 1, consistent with 
the fraction of sources expected to be below the SDSS 
limit (~ 40%) and therefore at high z. Only 33 of them 
have > 5<T detections in at least one PACS channel; these 
include the strongly lensed galaxies found by Negrello et 
al. (2010). 

5. LUMINOSITY FUNCTION OF HIGH-Z SUB-MM 
GALAXIES 

We have computed the LFs at rest-frame wavelengths 
of 100 and 250 /im in 4 redshift intervals (1.2 < z < 1.6, 
1.6 < Z < 2, 2 < z < 2.4, 2.4 < z < 4), exploiting the 
classical Schmidt's (1968) l/V maiX estimator, our redshift 
estimates, and if— corrections computed from the SEDs 
used for the redshift estimates. The V max was computed, 
for each galaxy, taking into account the redshift bound- 
aries of the bin and the maximum accessible redshifts 
implied by the 250 and 350 fjjn flux density limits. In 
the case of our preferred SED (SMM J2135-01012) the 
numbers of sources in these redshift bins are 900, 891, 
821, and 502, respectively. 



The estimated re dshifts 
http: / / people.sissa.it / ~gnuevo/photoz / 



available 



The error on z introduces both a statistical and a sys- 
tematic effect. The latter is related to the Eddington 
bias, that can be quantified with a Bayesian approach. 
The probability that the true redshift of a source is z 
when its estimated value is z e reads 



p(z\z e ) oc P(z e \z) n(z) 



(3) 



where P(z e \z) is the probability that the estimated red- 
shift of a source is z e when its true value is z and n(z) 
is the true redshift distribution. We take P(z e \z) to be 
a Gaussian with mean z and dispersion equal to the rms 
difference between spectroscopic and photometric red- 
shifts for the galaxies in Fig. [3] and n(z e ) as a proxy for 
n(z). The maximum likelihood estimate of z is the value 
for which p(z\z e ) peaks. All the luminosity functions 
presented in this paper are corrected for this bias. We 
caution, however, that this correction does not entirely 
account for the Eddington bias due to the errors on esti- 
mated luminosities. The effect is expected to be smaller 
than that of the error on z, which is the main source 
of the error on luminosity. A full treatment would re- 
quire extensive numerical simulations that are beyond 
the scope of this paper. 

As mentioned in §[2j we expect a residual contamina- 
tion of our reference sample by cold low-z galaxies with 
genuine r < 22.4 counterparts that could not be iden- 
tified with R > 0.8. To quantify the effect of this con- 
tamination on our luminosity function estimates we have 
exploited the H-ATLAS SDP sources with spectroscopic 
redshift in the Smith et al. (2011a) catalog. Out of a to- 
tal of 1096 galaxies with z spcc < 0.5, 682 comply with our 
selection criteria (Saso/an > 35mJy and S^o^m > 3<r). 
For these sources the use of the SMM J2135-0102 tem- 
plate for redshift estimates is inappropriate and yields a 
substantial positive offset and a rather broad dispersion 
of Zphot — Zspec (see §[4j) . In fact, after the Bayesian cor- 
rection described above, this template yields z p hot > 1-2 
for 160 of them. Of these, 88 fall in the first redshift bin, 
45 in the second, 21 in the third, and 6 in the fourth. 
If the sample of SDP galaxies with z spoc < 0.5 is rep- 
resentative and about 20% of the 5021 galaxies in our 
sample prior to the removal of R > 0.8 objects have 
unrecognized genuine SDSS counterparts, the low-z con- 
taminants would amount to 14% of sources in the first 
redshift bin, 7.4% in the second, 3.8% in the third, 1.8% 
in the fourth. On the other hand, it is estimated (Smith 
et al. 2011a) that w 5.8% of the 2418 R > 0.8 iden- 
tifications are spurious. If these spurious identifications 
have the same redshift distribution as the 3469 sources 
not identified with R > 0.8, we would be missing about 
4% of sources in each redshift bin. Finally, we have to 
take into account the incompleteness of the SDP cata- 
log. This can be done using the flux density dependent 
correction factors given in Table 2 of Rigby et al. (2011). 

The estimated LFs, taking into account all the correc- 
tions described above, are presented in Figs. [SHU and in 
Tables [1] and O The upper scale in these figures displays 
the SFR corresponding to the 100 or 250 /jm luminosity 
for the SMM J2135-0102 calibration giving 



L W0 , m _ 5Jxl0 23. SFR 



4 The results, however, are only weakly constrained by PACS 
data. We do not find significant differences if we use 5cr upper 
limits or ignore these data altogether. 



(4) 



WHz- 1 Moyr 1 ' 
while for £250^111 the coefficient is 1.4 x 10 23 . Since for 
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galaxies with intense star formation the rest-frame dust 
emission peaks in the range A « 90 — 100 /im, the 100 /im 
luminosity is a good estimator of the SFR. 

Note that the flattening of the LFs at the lowest lumi- 
nosities may be, at least in part, due to the overestimate 
of the accessible volume yielded by the 1/V max estimator 
for objects near to the flux limit, pointed out by Page & 
Carrera (2000; see also Eales et al. 2009). On the other 
hand, the Page & Carrera (2000) estimator holds under 
the assumption that the luminosity function varies little 
within the luminosity bin, while in our case it is very 
steep over most of the luminosity range. 

Fig. [B] illustrates the effect of varying the template 
SED, after correcting for the median offsets in the red- 
shift estimates highlighted by Fig. |3l The differences are 
quite limited. Such a stability of the LF estimates follows 
from some favorable circumstances. First, the large num- 
bers of sources in each redshift bin smooth out the effect 
of errors on redshift estimates. Second, for the redshift 
range considered here (1 < z < 4) one of the SPIRE 
wavelengths is always sampling directly the rest-frame 
SED in the range 100 < A < 125 jum, implying that K- 
corrections (and the related uncertainties) are minimal. 
Third, the LFs are only moderately sensitive to the un- 
certainty on redshift estimates. This can be shown as 
follows. The luminosity at a rest-frame frequency v is 
related to the flux density at the observed frequency z^ bs 
by 

L v = f Uohs jz ; k - , > (5) 



(l + z)xl ' 



where 



I 



L 



f obs + 



(6) 



Since, in our case, v w j/ bs (1 + z), I rs 1. Then it is 
easily checked that, in the redshift range of interest here, 
L v oc (1 + z) 3 , so that 



A log L 



Az 



ln(10) \l + z 



(7) 



The LF, <f>, is computed by weighting each galaxy by 
the inverse of the maximum accessible volume (Schmidt 
1968), i.e., 4> K At z > 1 we have, roughly, V oc 

(l + z) 



1 whence 



Alog< 



-1.2 



1 + 2 



(8) 



Since Az/(l + z) 0.2 the error on the LF determination 
due to uncertainty in the redshift estimate amounts to 
around 0.25 dex. The overall uncertainty is the sum in 
quadrature of the statistical error and of the error coming 
from the uncertainty on the redshift estimate. As shown 
in Fig. HI the correction for the Eddington bias is rather 
small compared to the overall uncertainty of the LF. 

In Figs. [5] and [7] our LF estimates are compared with 
those by Gruppioni et al. (2010) and Eales et al. (2010b), 
respectively. The latter authors adopted a grey body 
SED with T = 26 K and a dust emissivity index $ = 1.5. 
The assumed dust temperature is possibly more appro- 
priate for relatively low SFR/far-IR luminosity galaxies 
while typical dust temperatures of SPIRE-detected z w 2 



galaxies are somewhat higher (Chapman et al. 2010; Am- 
blard et al. 2010). The different choice for the SED has 
a substantial impact on the K— correction and therefore 
on the estimated rest-frame luminosity. For example, for 
a galaxy at z w 1.5 and a given observed 250 /im flux 
density, the Eales et al. SED yields a rest-frame 250 /im 
luminosity a factor « 1.5 higher than that obtained using 
either the SMM J2135-0102 or the Arp220 SED. Once we 
replace their if— correction with ours, the LF estimates 
by Eales et al. (2010b) are found to be in good agree- 
ment with ours (see Fig. [7|). The two estimates are to 
some extent complementary: the deeper samples used by 
Eales et al. have allowed them to reach lower luminosi- 
ties, while the larger H-ATLAS SDP area has allowed us 
to reach higher luminosities at high redshifts. 

Gruppioni et al. (2010) exploited the deep PACS 
data at 100 and 160 /im in the GOODS-N field (~ 
150 arcmin 2 ), obtained as part of the PACS Evolutionary 
Probe (PEP) survey, to estimate the 60 and 90 /im rest- 
frame LFs up to z ~ 3. They have used all the available 
data to derive the SEDs of their sources. 

About 31% of Eales et al. sources have spectroscopic 
redshifts and almost all the others have photometric red- 
shifts typically based on nine optical and near-IR bands. 
Gruppioni et al. (2010) have spectroscopic redshifts for 
~ 70% of their sources; for the remaining ~ 30% they 
have obtained photometric redshifts from multifrequency 
data. The agreement between our LFs and theirs is an 
additional confirmation that photometric redshifts de- 
rived from sub-mm photometry are good enough for the 
present purpose. 

In the sampled luminosity range, the LFs exhibit an ex- 
ponential fall off and a substantial luminosity evolution 
at least up to z w 2.5, while a weaker evolution at higher 
z is indicated (see Fig. |8]). The 250 /im luminosity corre- 
sponding to log(</>/Mpc~ 3 ) = -5 is log(L 2 5o/WHz _1 ) = 
26.02, 26.16, 26.32, and 26.38 for z « 1.4, 1.8, 2.2, 3.2, 
respectively. At 100 /im the corresponding luminosities 
are log^ioo/WHz" 1 ) = 26.65, 26.78, 26.94, and 27.01. 
These results are consistent with those based on PEP 
survey data (Gruppioni et al. 2010), which have poorer 
statistics at high z. We remark that our data supple- 
mented with those by Gruppioni et al. (2011) and Eales 
et al. (2010b) allow the determination of the LFs over 
only one order of magnitude in luminosity for 1 < z < 2, 
and over only a factor of about 3 for z > 2. 



6. 



CLUES ON STAR-FORMATION TIMESCALES IN 
MASSIVE HALOS 



We now discuss how the LFs of high redshift star 
forming galaxies at sub-mm wavelengths and the counts 
in sub-mm and mm bands concur with their clustering 
properties in probing the process of star formation in the 
progenitors of massive ETGs. 

6.1. Clustering properties and host halo masses of 
sub-mm galaxies 

Several lines of evidence indicate that high-z (sub-)mm 
bright galaxies are strongly clustered (Blain et al. 2004; 
Farrah et al. 2006; Magliocchetti et al. 2007; Viero et al. 
2009; Maddox et al. 2010; Hall et al. 2010; Cooray et al. 
2010; Dunkley et al. 2010; Amblard et al. 2011; Planck 
Collaboration 2011). The large-scale clustering power- 
spectrum probes the relation between the distribution 
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of visible sources and that of dark matter halos, which 
is a very sensitive function of the halo mass scale (e.g. 
Matarrese et al. 1997). 

The study of the angular correlation function of H- 
ATLAS SDP galaxies (Maddox et al. 2010) did not de- 
tect significant clustering for the 250 (im selected sample, 
while the clustering signal was found to be strong (al- 
though with large uncertainties) for samples selected at 
350 and 500 /im. The measurements are consistent with 
the idea that sub-mm sources consist of a low redshift 
population of moderately star-forming galaxies and a 
high redshift population of highly clustered star-bursting 
galaxies. The former, known to be weakly clustered (e.g. 
Madgwick et al. 2003; Guo et al. 2011; van Kampen et 
al. 2011), dominate the 250 /im selected sample, while 
the selection at longer wavelengths emphasizes strongly 
clustered galaxies with intense star-formation activity at 
z w 2 — 3, with typical halo masses » 1O 13 M , which 
are interpreted as the ancestors of present day massive 
ellipticals (Negrello et al. 2007). 

6.2. Timescales of sub-mm and UV bright phases 

The bright end of the LFs provides information on the 
average duty cycle of the star formation in the large halos 
indicated by the clustering analysis. This can be illus- 
trated using a simple model (hereafter referred to as the 
'toy model') relying on the following assumptions: 

• Most of the star formation occurs soon after the 
fast collapse phase of the halo, as identified by Zhao 
et al. (2003) and by several subsequent works (e.g., 
Diemand et al. 2007, Genel et al. 2010, Wang et 
al. 2011). 

• The halo formation rate is given by the positive 
term of the derivative with respect to the cosmic 
time of the Sheth & Tormen (1999, 2002) cos- 
mological mass function [eq. (1) of Lapi et al. 
(2006)]. This was shown to be a good approxi- 
mation at the redshifts (z > 1) and for the halo 
masses (Mh >3x 10 12 M ) of interest here (e.g., 
Haehnelt & Rees 1993; Sasaki 1994; see also the 
discussion by Lapi et al. 2006). 

• During the main episode of star formation, the SFR 
is proportional to the halo mass, i.e., SFR oc Mjj, 
with a small scatter. 

• The duration At s f of the main star formation 
episode, before quenching by the AGN feedback, is 
roughly constant for the considered range of halo 
masses (3 x 10 12 M < M H < 3 x 10 13 Af Q ) and 
shorter than the expansion timescale at the source 
redshift. 

This very simple model provides a good fit to the data 
(Figs. [6]and El) if At s{ w 7 x 10 8 yr and: 



SFR = 3£ 



M, 



H 



10 12 M, 



1 



2.5 



M yr" 



(9) 



with a rather small (±35%) scatter. We have used Eq. @ 
with the calibration appropriate for SMM J2135-0102 to 
pass from star formation rate to FIR luminosity, and the 
SED of this object to go from the FIR luminosity to 



monochromatic luminosities. The fit would significantly 
worsen if the halo mass range associated to individual 
galaxies extends above M# > 3 x 1O 13 M . This sug- 
gests that the star formation becomes very inefficient in 
the most massive halos, possibly due to the cooling time 
becoming longer than the expansion time. The bright 
end of the LFs is quite sensitive to the value of erg; e.g., 
raising erg from our fiducial value of 0.81 to 0.9 would 
increase by a factor m 2 the 250 ^m LF at z « 2.5 and 
~ 3x 10 26 W Hz -1 . Interestingly, the duration of 
the main star-formation episode yielded by the toy model 
matches that inferred from the observed a-enhancement 
of massive local ETGs (see §[TJ). 

A rough estimate of the mass in stars at the end of 
the main star-formation episode, for given halo mass and 
redshift, can be obtained multiplying the appropriate 
SFR by its duration Ai s f « 7 x 10 s yr, and correcting 
for the fraction of stellar mass returned to the interstel- 
lar medium. For a Chabrier IMF the latter amounts to 
~ 30% after 0.1 Gyr from a burst of star formation, and 
increases to « 40% after 1 Gyr and to 50% after several 
Gyrs. Combining with the halo mass function at that 
redshift we get an estimate of the stellar mass function 
of ETGs. A comparison with observational determina- 
tions at different redshifts is shown in Fig. [51 The agree- 
ment is remarkably good, considering the crudeness of 
the approach. This suggests that the toy model captures 
the basic aspects of the star formation in massive galaxy 
halos at high redshift. 

A comparison of our LFs with the UV ones at z w 2 — 3 
(Sawicki & Thompson 2005; Reddy & Steidel 2009) for 
galaxies with comparable SFRs shows that the UV space 
densities are a factor > 100 lower, implying that the 
UV bright, dust-free phase is much shorter than the far- 
IR bright phase. The far steeper slope of the UV LF, 
compared with the sub-mm ones, for SFR ~ hundreds 
M yr _1 implies that, for the corresponding range of 
galaxy masses, the UV-bright phase is shorter for larger 
SFRs, i.e. for more massive galaxies, as predicted by 
Mao et al. (2007). 

If, as implied by most current semi-analytical models 
(Granato et al. 2001, 2004; Croton et al. 2006; Hopkins 
et al. 2008; Somerville et al. 2008), the black hole growth 
is linked to star- formation, most of the present-day black 
hole mass should have been accreted by the end of the 
star-forming phase (see Bonfield et al. 2011). Therefore, 
the local ratio M./M* « 2 x 10~ 3 (Marconi & Hunt 2003) 
between the black hole mass M. and the stellar mass M* 
of the host ETG should be already in place then (but see 
also McLure et al. 2006 and Peng et al. 2006). If so, the 
stellar mass function translates immediately into a black 
hole mass function which, in turn, can be translated into 
a bolometric luminosity function assuming that the black 
holes, immediately before being switched off, emit at the 
Eddington limit for a time i v i s - We get: 



L 



BH,bol 



79 



Mr, 



(10) 



The corresponding B— band luminosity is Lg — 
0.1(10//cboi)iBH,boi, where fcboi is the bolometric correc- 
tion (Marconi et al. 2004; Hopkins et al. 2007). We 
can then compare the comoving density of galaxies with 
that of the associated QSOs. For example, at z « 2 
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there are * 2.7 x 10 5 sub-mm galaxies per Mpc 3 with 
M* > 2 x 10 11 M , while the density of QSOs brighter 
than the corresponding luminosity of Lb ~ 1.6 x 10 12 Lq 
(fcboi « 10) is » 8 x 10- 7 Mpc -3 (Croom et al. 2004). 
This implies that the optical visibility time of the QSOs 
is about a factor of 30 shorter than the duration of the 
sub-mm bright phase, in agreement with independent es- 
timates (Shankar et al. 2004; Marconi et al. 2004; see 
also Lapi et al. 2006). 

The physical model by Granato et al. (2001, 2004), 
further elaborated by Lapi et al. (2006), complies with 
the indications coming from the previous analysis. In 
this model the star formation in massive galaxies occurs 
on a timescale of several 10 8 yr, is very soon obscured by 
dust and is stopped by quasar feedback. The star forma- 
tion is triggered by the rapid cooling of the gas within 
a region with an approximate size w 1/3 — 1/4 of the 
halo virial radius, i.e., 20 — 30 kpc, encompassing about 
40% of the total mass (dark matter plus baryons), and 
is regulated by the energy feedback from SNe and AGNs 
(the latter being relevant especially in the most massive 
galaxies). As a result, the star formation is a very ineffi- 
cient process, as proved by the fact that on the average 
only about 5 — 10% of the baryons in the Universe are 
converted into stars (Fukugita & Peebles 2004), with a 
possible maximum of 15—20% for halo masses ~ 10 12 Mq 
(Shankar et al. 2006; Moster et al. 2010). It is worth 
noticing that the Negrello et al. (2007) predictions for 
the lensed sub-mm bright galaxies follow directly from 
the LF and redshift distribution yielded by the Lapi et al. 

(2006) model. We have updated this model by replacing 
the Arp220 GRASIL SED used, e.g., by Negrello et al. 

(2007) with the SED of SMM J2135-0102. The resulting 
LFs are shown in Figs. [5] and [71 labeled as 'Full model'. 
We stress that the LFs from the Full model are not fits 
to the data, but have been computed with the same pa- 
rameters used by Lapi et al. (2006; see their Table 1). 
This model has proven to be successful in reproducing 
a wealth of observables, including high redshift quasar 
LFs. 

Note that, in the present framework, sources dominat- 
ing the high-z LF are proto-spheroidal galaxies that are 
in passive evolution at z < 1.5 — 1 and therefore do not 
contribute to the low-z FIR/sub-mm LFs computed, e.g., 
by Dye et al. (2010), Vaccari et al. (2010), and Dunne 
et al. (2011). 

An alternative model whereby intense star-formation 
at z rs 2 — 3 is not supported by mergers was proposed 
by Dekel et al. (2009). In this model, star-formation is 
driven by steady cold streams supplying gas at an ap- 
proximate rate of 

(Hi 

where /o.i65 is the baryonic fraction in the halos in units 
of the cosmological value, /& = 0.165. The dependence of 
the cold gas inflow rate, M, on Mh and z are similar to 
those implied by our toy model [Eq. ©], implying that 
the SFR must be roughly proportional to M. To account 
for the SFR indicated by the Herschel data, the star- 
formation efficiency must be very high, ~ 70% for Mh » 
10 12 Mq and z ~ 1.5. For comparison, the mean cosmic 



star-formation efficiency is ^star/^baxyon ~ (5-8 ± 1.1)% 
(Fukugita & Peebles 2004; Komatsu et al. 2011). In 
halos more massive than w 10 13 M© cold streams are 
suppressed by shock heating. 

7. (SUB)-MM COUNTS AND REDSHIFT DISTRIBUTIONS 

Figs. [TU1 QTJ and [T^] compare the predictions of the 
'full model' with the observed counts in Herschei/SPIRE 
bands. In these figures the contributions to the counts 
of proto-spheroidal galaxies have been complemented 
with those of normal late-type and starburst galaxies 
computed by Negrello et al. (2007). As for mas- 
sive proto-spheroidal galaxies, the main difference with 
the counts in the latter paper comes from having re- 
placed the Arp220 SED yielded by GRASIL with that of 
SMM J2135-0102. The contribution from lensed proto- 
spheroidal galaxies has been estimated using the amplifi- 
cation distribution of Perrotta et al. (2003) and Negrello 
et al. (2007). 

An interesting prediction of the model is that massive 
proto-spheroidal galaxies dominate the (sub-)mm counts 
over a limited flux density range (about a decade). At 
250 fim the Euclidean normalized differential counts of 
these objects peak at » 30 mJy, i.e. roughly at the 
detection limit of the H-ATLAS survey, and sink down 
rapidly at fainter fluxes where the contribution of star- 
burst galaxies takes over and accounts for the results of 
the P(D) analysis by Glenn et al. (2010). Since the 
proto-spheroidal galaxies are mostly at z > 1.5 while 
late-type/starburst galaxies are mostly at z < 1.5, the 
model implies that the redshift distribution drifts to- 
wards lower redshifts as we go fainter. This is not in 
contradiction with the finding by Bourne et al. (2011) 
that the contribution of z < 0.28 galaxies to the cosmic 
infrared background is < 5% in the SPIRE bands. 

The model accurately fits the source counts from 
250 /zm to 850 /im (Figs.[IDl-fI5|. The observed counts at 
1.1 mm (Fig.rjJl span the peak of the Euclidean normal- 
ized counts, while those at 1.4 and 2 mm fFigs.[T5land [T6|) 
only cover the brightest tail of the counts, dominated by 
strongly lensed sources (apart from radio sources). The 
model somewhat overestimates the counts of strongly 
lensed galaxies at 1.4 and 2 mm, and the most recent 
counts at 1.1 mm. This may suggest that higher- z galax- 
ies, that yield larger and larger contributions to the 
bright counts at increasing mm wavelengths, have SEDs 
slightly colder than SMM J2135-0102 and closer to that 
of G15.141. 

The flux density range over which massive proto- 
spheroidal galaxies dominate the counts increases some- 
what with increasing wavelength (Figs. fT0l- fT6)) . reflect- 
ing the increase of the redshift range through which the 
strongly negative K— correction makes the flux corre- 
sponding to a given luminosity essentially independent 
of distance. The broadening of the peak in the Euclidean 
normalized counts of proto-spheroidal galaxies is not very 
large, however, because massive halos become increas- 
ingly rare at high z. At all wavelengths the bright counts 
of proto-spheroidal galaxies drop down very steeply, re- 
flecting the exponential decline of the halo mass function. 
This rapid fall-off, borne out by the data, can hardly be 
accounted for by phenomenological models that evolve 
the LF of local populations of dusty galaxies backwards 
in time: spheroidal galaxies are essentially in passive evo- 
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lution since z « 1.5 and are therefore not represented in 
local LFs at far-IR to mm wavelengths. 

The physical model of Lacey et al. (2010) works pretty 
well at 850 /im but does not correctly predict the rapid 
fall-off of the counts above « 30 mJy at 250 /im, 350 /im 
and 500 /jm (Figs. IT0l-fT3|. The decline of the counts is 
most easily understood if the duration of the most active 
star- formation phase in massive halos is relatively short, 
as also required by the observed a-enhancement (see §[J). 
However, in merger-driven evolutionary models, such as 
the one by Lacey et al. (2010), the star formation does 
not truncate after < 1 Gyr. At 850 /jm consistency with 
the data can be recovered if sub-mm bright galaxies con- 
tain large amounts of cold dust, but this leads to 250 fxm, 
350 /im and 500 fim counts less steep than observed. 

The fact that using, for high-z galaxies, the SMM 
J2135-0102 SED we are able to simultaneously reproduce 
the source counts over a broad wavelength range indi- 
cates that this SED is indeed typical for sub-mm bright 
galaxies. In fact the results on the counts vary substan- 
tially if we use different SEDs (see Fig. [17]). Almost by 
construction, in all cases the 250 /im (the main selection 
wavelength) counts are accurately reproduced. But if, 
for example, we use the SED of G15.141 instead of the 
SMM J2135-0102 one, the longer wavelength counts are 
somewhat underestimated. The situation is better, but 
still not as good as for our preferred SED, if we use the 
Arp220 SED. The quite cold (T d = 26 K) SED used by 
Eales et al. (2010b) leads to counts substantially in ex- 
cess of the observed ones. Thus the 'colors' of the counts 
constrain the SEDs of high redshift star forming and in 
particular on the abundance of galaxies with cold dust 
temperatures. As pointed out by Hwang et al. (2010) 
apparently very cold values of Td could be caused by an 
overestimation of the sub-millimeter fluxes due to blend- 
ing problems or inappropriate single-temperature SED 
fitting. 

Fig. HU compares the redshift distributions of z > 1 H- 
ATLAS SDP sources detected at S > 35, 40, 48 mJy at 
250, 350 and 500 /im, respectively, with those predicted 
by the model. The effective flux correction factors tabu- 
lated by Clements et al. (2010) have been applied for this 
comparison. Since the redshift distributions are a key in- 
gredient for the successful prediction of the abundance of 
strongly lensed galaxies and of their redshift range (Ne- 
grello et al. 2010), they cannot be badly wrong. Thus, 
the agreement between our estimates and model predic- 
tions is a confirmation of the global consistency of our 
results. 

The shift in the redshift distribution of sub-mm 
bright galaxies to higher redshifts with increasing selec- 
tion wavelength happens because the strongly negative 
K— correction extends to larger and larger redshifts and 
the lensing probability increases with z. For example, we 
expect that, after removing local objects that are found 
to be only k 1/3 of the total (Vieira et al. 2010), the 
redshift distribution of sources with Si. 4 > 10 mJy has a 
broad peak around z 4.5. Using the SMM J2135-0102 
SED we also find that the observed 1.4 — 2 mm spectral 
indices are correlated with z: for z < 4.5, af 4 > 3.2, 
while for z > 4.5, a\ 4 < 2.6. Thus the value of this 
spectral index may be a rough redshift indicator. 

In Fig. [19] we illustrate the contribution of proto- 



spheroids to the extragalactic infrared background in- 
tensity (Lagache et al. 1999; Hauser & Dwek 2001) 
according to our physical model. At A > 850 /im the 
background is almost entirely accounted for by high- 
redshift proto-spheroids, while these contribute about 
65% and about 50% of the background at A « 500 /^m 
and A ~ 250 /im, respectively. 

8. CONCLUSIONS 

We have exploited the H-ATLAS SDP survey data to 
investigate the evolution of the 100 and 250 /am LFs of 
bright star-forming galaxies (SFR > 10 2 M Q yr~ 1 ) at 
z > 1. Redshifts have been estimated using 3 SED tem- 
plates representative of the range of well measured SEDs 
of galaxies with such a high star formation. The rms 
uncertainties on redshift estimates have been assessed 
comparing our redshift estimates with spectroscopic red- 
shift measurements for 39 H-ATLAS galaxies at z > 0.5 
as well as by means of simulations. Both methods yield 
rms values ofAz/(l + z)«0.2or smaller. The LFs de- 
rived using the redshift estimates based on each of the 
3 SED templates are very similar to each other. The 
uncertainties due to the spread of redshift estimates are 
added in quadrature to Poisson errors to compute the 
global uncertainties on the LFs. Our LF estimates are 
in close agreement (in the common redshift and luminos- 
ity range, after applying the same if —corrections) with 
those at 250 /im by Eales et al. (2010b), based on a sub- 
stantial fraction of spectroscopic redshifts complemented 
with photometric redshifts of optical/near- IR counter- 
parts, as well with those at 90 /im by Gruppioni et al. 
(2010), based on an even higher fraction of spectroscopic 
redshifts. This is a further confirmation that our red- 
shift estimates are sufficiently accurate to allow reliable 
estimates of the LFs. 

The SED of SMM J2135-0102 was found to perform 
significantly better than the others in the tests we made, 
at least up to z 3, and was therefore adopted as our 
reference. Remarkably, this SED allowed us to simulta- 
neously fit the counts over a broad wavelength range. 

We find (see Fig. [5]) a significant luminosity evolution 
at least up to z « 2.5 while the luminosity function 
shows a modest variation between the last two redshift 
bins, centered at z « 2.2 and z « 3.2, respectively, even 
though the corresponding time interval is « 1 Gyr, sub- 
stantially larger than the time interval (« 0.65 Gyr) be- 
tween the central redshifts of the second and third bin 
(z « 1.8 and w 2.2). This is consistent with the results 
based on PEP survey data (Gruppioni et al. 2010). We 
show that the evolution of the LF reflects that of the 
halo formation rate if, for the very massive galaxies rep- 
resented in our sample, the star-formation rate obeys a 
simple relationship with halo mass (estimated from the 
clustering properties) and redshift [Eq. (|9])] and the life- 
time of the main star-formation phase is Af s f « 7 x 10 s 
yr, consistent with the constraint coming from the ot- 
cnhanccmcnt observed in the most massive ETGs. 

The stellar mass function resulting from the derived 
star-formation rates and timescale At s f nicely fits the ob- 
served stellar mass function of passively evolving ETGs 
at-rwl-2. 

A comparison of sub-mm and UV LFs indicates that 
the UV visibility time of massive galaxies is much (> 100 
times) shorter than the duration of the sub-mm bright 
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phase (At s { «7x 10 s yr). This implies that dust forms 
very rapidly after the onset of the main episode of star 
formation, either in the ISM (Draine et al. 2009; Dunne 
et al. 2011) or by effect of Type-II SNae (Dunne et al. 
2003, 2009; Matsuura et al. 2011). 

As discussed by Mao et al. (2007), a longer UV bright 
phase is expected for less massive galaxies. This predic- 
tion is supported by the flatter slope of the far-IR LF, 
compared to the UV one, for SFR ~ hundreds M Q yr -1 , 
implying a decrease with increasing galaxy mass of the 
ratio of UV to far-IR space densities, i.e. of the ratio 
between UV and far-IR lifetimes. 

In the same vein we find that the duration of the op- 
tically bright QSO phase is of order 1/30 of At s {, i.e., of 
about 2 — 3 x 10 7 yr, consistent with independent evi- 
dences. On the other hand, the bright end of the QSO 
LF is somewhat less steep than that of sub-mm galaxies. 
This is not surprising because, although the evolution of 
the SFR and of accretion into the central black hole may 
well be linked, the relationship is mediated by several 
steps that introduce a substantial scatter that, coupled 
with the curvature of the LF, flattens its bright tail. 

The 100 and 250 /xm LFs at different redshifts are quite 
well reproduced by the physical model of ETG formation 
and evolution by Granato et al. (2001, 2004), further 
elaborated by Lapi et al. (2006), without any adjust- 
ment of the parameters. As discussed in these papers, 
the model is built in the framework of the standard hi- 
erarchical clustering scenario. Many simulations (e.g., 
Zhao et al. 2003, Diemand et al. 2007, Genel et al. 
2010, Wang et al. 2011) have shown that the growth 
of a halo occurs in two different phases: a first regime 
of fast accretion in which the potential well is built up 
by the sudden mergers of many clumps with comparable 
masses; and a second regime of slow accretion in which 
mass is added in the outskirts of the halo, only occa- 
sionally affecting the central region where the galactic 
structure resides. According to the model, the fast ac- 
cretion phase triggers a burst of star formation that, in 
massive halos at z > 1, starts an evolutionary sequence 
that can be summarized as follows. There is an early, 
short phase of (almost) dust-free star formation when the 
galaxy shines as a bright UV source. It is followed by a 
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dust-enshrouded star formation phase when the galaxy 
shines in the far-IR/sub-mm range. The duration of both 
the UV bright phase and the far-IR/sub-mm bright phase 
is shorter for the most massive galaxies, with the highest 
SFRs; for these objects the UV phase lasts < 10 7 yr and 
the far-IR/sub-mm phase lasts < 10 9 yr. Then there is 
a phase, lasting several 10 7 yr, when the nucleus shines 
as a bright QSO after having swept away most of the in- 
terstellar gas and dust. Finally, passive evolution of the 
stellar populations follows, and the galaxy evolves into a 
local ETG. 

According to this model, star-forming proto-spheroidal 
galaxies account for a substantial fraction of the cosmic 
infrared background (see Fig. and dominate the cos- 
mic SFR at z > 1.5, while at lower z the SFR is domi- 
nated by late-type (normal and starburst) galaxies. This 
model was the basis for the successful predictions the 
sub-mm counts of strongly lensed galaxies by Perrotta et 
al. (2003) and Negrello et al. (2007). It also accurately 
reproduced the epoch-dependent galaxy luminosity func- 
tions in different spectral bands, as well as a variety of 
relationships among photometric, dynamical and chemi- 
cal properties, as shown in previous papers (see Table 2 
of Lapi et al. 2006 and additional results, especially on 
the galaxy chemical evolution, in Mao et al. 2007 and 
Lapi et al. 2008). 
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Figure 1. Optical (SDSS r— band) magnitudes as a function of rcdshift for several SEDs, normalized to Sssuum = 33.5 mjy, the 5<r 
detection limit for H-ATLAS SDP galaxies. The red curves refer to the mean SEDs of optically identified z < 0.5 SDP galaxies in the 
lowest [9.5 < log(L(just/ Lq) < 10; solid line] and in the highest [11.5 < log(L^ ust / Lq) < 12; dashed line] luminosity bins of Smith ct al. 
(2011b). The blue curves refer to the same sample by Smith et al. (2011b) but for the lowest [—11.5 < log(sSFR) < —11; dashed] and for 
the highest [—8.5 < log(sSFR) < —8; solid] specific star formation rate (in yr — x ) bins. The dashed black line refers to the SED of Arp220 
(a local ULIRG) and the black solid line to that of SMM J2135-0102 ('The Cosmic Eyelash' atzss 2.3). As explained in § 2, the figure is 
meant to illustrate that the high-2 H-ATLAS sources must have global SEDs different from the low-2 ones. The latter are relatively cool 
and unobscured, while the high-2 population is dominated by more obscured objects with higher star-formation rates and warmer dust 
temperatures. 



14 



A. Lapi et al. 



SEDs 



i i i i i i i i I i i i i i i i i i I i i i i i i i i i I i i i i i i i i i I i i i i i i i i i I i i i i i i i i i 



E 
=1 

o 
o 

@ 



cn 
O 



• SMM J2135-0102 

lvison+ 1 
Swinbank+ 1 




■5E 
■6 



SMM J2135-0102 (GRASIL) 

-■- Arp 220 (GRASIL) 
- - G1 5.141 (2GB) 

I I I I I ' 



1 







1 2 3 

log A [/i.m] 



4 



Figure 2. The SED of SMM J2135-0102 ('The Cosmic Eyelash') as modeled by GRASIL (blue solid line); data are from Ivison et al. 
(2010) and Swinbank et al. (2010). The SED of Arp220 (red dot-dashed line) as modeled by GRASIL, and that of G15.141 modeled as a 
double-temperature grey-body (green dashed line; see text for details) are also shown for comparison. All SEDs are in the rest-frame and 
are normalized at 100 (mi. 
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Figure 3. Photometric rcdshift estimates based on the SEDs of (from top to bottom) SMM J2135-0102, Arp220, G15.141, and the 
average low-z SED from Smith et al. (2011) compared in terms of the quantity Az/(1 + z) = (zp^ot — z S pec)/(l + z spec ) with spectroscopic 
measurements for the strongly lensed galaxies of Negrello et al. (2010; red filled circles) and Cox et al. (2011; purple filled square), and 
for the H-ATLAS SDP galaxies with z > 0.5 of Bonfield et al. (2011; blue asterisks) and Smith et al. (2010; black dots). The dashed lines 
correspond to a 20% difference. 
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Figure 4. Comparison of the redshift distributions of sources in our reference sample (S250/jm > 35mjy, S^o^m > 3<r, and no R > 0.8 
optical identifications; 3469 sources, i.e. m 50% of the 6876 H-ATLAS SDP galaxies; the selection has mostly excluded the z < 0.5 sources, 
and in fact there is another peak in the redshift distribution of the total SDP sample at z < 0.5, see Smith et al. 2011a) estimated using 
the SEDs of Arp 220 (green dashed line), G15.141 (blue dot-dashed line), and of SMM J2135-0102 (black solid line), correcting for the 
offsets estimated from Fig. [3] For our preferred SED, the one of SMM J2135-0102, 3093 are estimated to be at z > 1. Correcting for the 
contamination by cold low-z galaxies with genuine r < 22.4 counterparts not identified with R > 0.8 (see §[5jl, we estimate that the fraction 
of SDP galaxies with S250;jm > 35mjy (6100 in total) that arc at 2 > 1 is of ?s 45%. 
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Figure 5. Luminosity functions at 100 (im in four redshift intervals compared with models described in the text (solid line: full model; 
dashed line: toy model). The flattening of the LFs at the lowest luminosities may be, at least in part, due to the overestimate of the 
accessible volume yielded by the 1/V m ax estimator for objects near to the flux limit (Page & Carrera 2000). The 90 fim LFs (crosses) 
estimated by Gruppioni et al. (2010) using PACS data are also shown for comparison. 
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Figure 6. Comparison of the luminosity functions at 100 fim computed using different SED templates - SMM J2135-0102 (blue circles), 
Arp220 (red squares), and G15.141 (green stars) — after correcting for the median offsets in the redshift estimates highlighted by Fig. [3] 
To ease the comparison, only statistical uncertainties are reported. 
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Figure 7. Luminosity functions at 250 (im in four redshift intervals based on our reference sample, compared with the predictions of the 
models discussed in the text (solid line: full model; dashed line: toy model) and with the estimates by Eales et al. (2010b), after applying 
the A"— corrections based on the SMM J2135-0102 SED, for common redshift ranges. 
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Figure 8. Synoptic view of the 100 fira and 250 fim LFs for our reference sample in four redshift intervals. In the 100 fim panel, open and 
filled symbols show our estimates before and after the correction for the Eddington bias, respectively (see § 4 for details). 
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Figure 9. Stellar mass functions of ETGs at different redshift predicted by the toy model compared with observational determinations. 
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Figure 10. Comparison of the observed Euclidean normalized counts at 250 fim (Clements et al. 2010; Oliver et al. 2010; Glenn et al. 2010; 
Bethermin et al. 2010) with the predictions of our full model (solid line) and of the semi-analytic model of Lacey et al. (2010; crosses). The 
different lines show the contributions to the counts of the various populations included in our full model: unlensed (blue dot-dashed line) and 
strongly lensed (purple dotted line) proto-spheroidal galaxies (host halo masses are within the range 2.5 X 10 11 M < M H < 3 X 10 13 M ), 
normal late- type plus starburst galaxies (green dashed line), radio sources (magenta triple-dot-dashed line). The radio source counts are 
from the De Zotti et al. (2005) model. The grey shaded area shows the counts of local galaxies estimated by Serjeant & Harrison (2005; 
SH05). 
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Figure 11. Same as in Fig. HOI but at 350 fim. 
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Figure 12. Same as in Fig, HOI but at 500 fim. Note that the flux densities in the Table 1 of Clements et al. (2010) must be multiplied by 
the correction factors given in the same table to produce the correct flux densities; the lowest corrected 500 micron flux density is 40 mjy. 
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Figure 13. Comparison of the observed integral counts at 850 (im (Zemcov et al. 2010; Knudsen et al. 2008; Coppin et al. 2006) with the 
predictio ns o f our full model (solid line) and of the semi-analytic model of Lacey et al. (2010; crosses). The lines have the same meaning 
as in Fig. 1101 
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Figure 14. Comparison of the observed Euclidean normalized counts at 1.1 mm (Scott et al. 2010; Perera et al. 2008; Austermann et al. 
2009, 2010; Hatsukade et al. 2011) with the predictions of our full model (solid line). The lines have the same meaning as in Fig. 1101 
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Figure 15. Comparison of the observed Euclidean nor mali zed counts at 1.4 mm (Vieira et al. 2010) with the predictions of our full model 
(solid line). The lines have the same meaning as in Fig. 1101 
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Figure 16. Comparison of the observed Euclidean nor mal ized counts at 2 mm (Vieira et al. 2010) with the predictions of our full model 
(solid line). The lines have the same meaning as in Fig. 1101 
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Figure 17. Constraints on the SED of high-z galaxies from multi-wavelength source counts. We have determined the luminosity functions 
at different redshifts > 1 using our reference sample and the SEDs of Arp220 and G15.141, in addition to SMM J2135-0102, and used 
them to compute the number counts at several wavelengths, adding the contributions of (lower-z) spiral and starburst galaxies as given by 
the Negrello et al. (2007) model. In all cases we recover a good fit at 250/^m, essentially by construction since this is the main selection 
wavelength. At longer wavelengths, however, while the counts obtained using the SMM J2135-0102 SED are still nicely consistent with the 
data, the other SEDs do not yield good fits. The counts given by the luminosity functions estimated by Eales et al. (2010), K-corrected 
with their SED, are also shown for comparison. In this case too the observed 250/xm counts are accurately reproduced. 



30 



A. Lapi et al. 



Submm red s h if t distribution 



S>35 mjy @ 250 fim 




1 2 3 4 5 6 



S>40 mjy @ 350 

4 I 




1 2 3 4 5 6 



S>48 (40) mjy ® 500 fjm 




1 2 3 4 5 6 



Figure 18. Comparison between the estimated redshift distributions (histograms) of proto-spheroidal galaxies at 250/^m (SuBOttm. > 
35 mjy), 350^tm (S350fim > 40 mjy), and 500 (im (SsoO/jm > 48 mjy), and the predictions of our full model, using the same parameters 
as in Lapi et al. (2006); solid lines refer to unlensed, while dashed lines also include the contribution of lensed, proto-spheroids. Because 
of the flux boosting due to confusion, the catalogued 500 fim flux densities have to be scaled down by the factors given by Clements et 
al. (2010). For example, the 500 fim flux density limit of 48 mjy for catalogued sources corresponds to a boosting corrected flux limit of 
40 mjy. As discussed in the text (see also the caption to Fig. |4}, we estimate that pa 55% of SDP galaxies with S^sOum > 35 mjy are at 
z < 1. Such galaxies belong to a different population comprising late-type normal and starburst galaxies. 
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Figure 19. Contributions of high-z proto-spheroidal galaxies to the (sub-)mm extragalactic background, according to our model, compared 
with the observational estimates of the total background intensity (Lagache et al. 1999). ETGs account almost entirely for the background 
at A > 850 (im. Contributions from lower redshift (z < 1 — 1.5) late type, normal and starburst galaxies become increasingly important at 
shorter wavelengths. 
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Note. — Based on our reference sample (S250fim > 35 mjy, 535o^ m > 3<x, and no R > 0.8 
optical identifications). Rcdshift estimates relying on the SMM J2135-0102 SED. Only statistical 
uncertainties are quoted; an additional uncertainty of A log (f> xs 0.25 related to the rcdshift estimate 
has to be summed in quadrature. 
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